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Surface plasmon-polaritons in a multilayer structure consisting of a metallic film and one or more layers of 
nanocrystal (NC) quantum dots (QDs) are studied. It is shown that there is resonance coupling between the 
plasmon-polaritons propagating along the metal/NC-layer interface and excitons confined in the dots, which 
makes a considerable effect on the optical properties of the structure unless the dispersion of the QD size is too 
large. This coupling can be explored in order to selectively excite QDs of different size by using an attenuated 
total reflection (ATR) structure. It opens the possibility of control of the relative intensity of light of different 
color, emitted by the QDs of different size. 
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INTRODUCTION 

Colloidal chemistry techniques have proved the ability to 
sinthesize high quality nanocrystals (NCs) of II- VI semicon- 
ductors, possesing the properties of quantum dots (QDs). 
With these NCs, a broad variety of nanostructures can be 
prepared, including multilayer structures of QDs of different 
average size, deposited on different substrates |1|. Combin- 
ing such structures with other materials, such as organic di- 
electrics 12 or graphene |3|, can result in new interesting 
physics and applications related to the energy transfer from 
photoexcited NCs to these materials. Recent experiments |4| 
demostrated the existence of a strong coupling between exci- 
tons confined in nanocrystal quantum dots (QDs) and surface 
plasmon-polaritons (SPPs) propagating along the interface of 
a silver film and the QD layer deposited on top of it. The ob- 
servation was achieved by measuring attenuated total reflec- 
tion (ATR) spectra of the structure. In this work we develop 
the theory of the observed effect taking into account the the 
QD size dispersion. We will show that the resonance SPP 
coupling excitons confined in QDs can be considerable unless 
the dispersion of the QD size is too large or the dots are too far 
from the metal/dielectric interface. The resoanance coupling 
can be used for controllable pumping the dots in order to ex- 
plore their uncomparable luminescence properties |[l]. Based 
on the calculated results, we will discuss possible applications 
of ATR structures containing, beyond the silver film, several 
layers of NCs of different average size. 



DIELECTRIC FUNCTION OF A COMPOSITE MATERIAL 
CONTAINING QUANTUM DOTS 

A single QD can be described by the electronic susceptibil- 
ity taking into account inter-band transitions: 
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where the sum runs over confined exciton states with ener- 
gies Ei (i = denotes exciton vacuum), c?cu is the transition 



dipole moment matrix element between valence and conduc- 
tion bands of the underlying semiconductor material, V is the 
QD volume and 
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with gi and (fe, r^) denoting the degeneracy factor and 

the wavefunction of the corresponding exciton states. 
The quantity 
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can be regarded as a QD dielectric function, where is the 
high-frequency dielectric constant of the QD material. How- 
ever, it can be more convenient to define the QD polaiizability 
(assuming its spherical shape), 
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where a is the QD radius and is the dielectric constant of 
the host material. Using these quantities (2j3i, one can cal- 
culate an effective dielectric function, e*, of the composite 
material containing uniform size QDs embedded in the host 
matrix using one of the following schemes: 

(i) Maxwell-Garnett approximation (MGA) valid in the low 
QD concentration limit 0, 
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(ii) Bruggemann mean field approximation (BA) lH, 
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where / — ^Na^ is the volume fraction of QDs. BA is 
valid when / ^ 0.5. MGA can be extended up to / ^ 0.1 
using a renormalized polarizability which takes into account 
dipole-dipole interactions between QDs as explained in Ref. 
7 

The modified MGA (MMGA) proposed in in Refs. El [8] 
allows for taking into account the QD radius dispersion. Let 
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F{a) denote the radius distribution function. The renormal- 
ized polarizabiHty, a*, can be calculated by the following 
equations, 

= / ^(-')^-^^^^J^Wda' ; (6) 



e(a') = ^iVa(a') / ^^^^F{a)da (7) 
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where a{a) is the polarizabiHty of a QD of radius a. For 
uniform-size QDs Eqs. (j6]) and (j?]) reduce to 



fa 

where a — a/a^. a* replaces a in Eq. (j4]l [notice that for 
small QD volume fraction fa* ^ a]. 



THE CASE OF CDSE SPHERICAL QDS 

Within the simplest model neglecting the Coulomb interac- 
tion between the electron and hole (strong confinement limit) 
and multiple sub-band structure of the valence band, the QD 
exciton spectrum is given by |9], 
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where Eg is the band gap energy (=1.75 eV for CdSe), /i is 
the electron-hole reduced mass and = tt, ^2 ~ 4.49, ^3 « 
5.76, ... are the roots of the spherical Bessel functions. At the 
same time, /3i = 2 for all i > 1. The dipole moment matrix 
element can be expressed as 
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where mp is the free electron mass and 2p1^/mQ ^20 eV. 
We take €h — 1.5 for PMMA matrix and assume a Gaussian 
distribution of QD radius. 
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with a and Aa being the average radius and the dispersion, 
respectively. The real and imaginary parts of the effective di- 
electric function of the QD/PMMA composite calculated for 
/ = 0.1 is shown in Fig. 1. A small imaginary part, F = 1 
meV was introduced in the denominators of Eq. ([TJ. 

As it can be seen from Fig. 1, the imaginary part of the 
effective dielectric function has a resonance at the exciton fre- 
quency, with the broadening depending on the QD concentra- 
tion and size dispersion. The real part of e* can reach negative 
values in the vicinity of the exciton frequency only for vanish- 
ing dispersion [Fig. [TJ^)]. while for larger dispersion it always 
remains positive [compare with Figs. [TJb) and[T|c)]. 
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no. 1: (color on line) Real (solid lines) and imaginary (dashed 
lines) parts of the CdSe QD - PMMA composite dielectric function 
for average QD radius a = 3nm and filling factor / = 0.1, calcu- 
lated using MMGA with QD radius dispersion Aa = (panel a), 
Aa = 0.005a (panel b), and Aa = 0.1a (panel c). 



SPP AT THE METAL/NC-LAYER INTERFACE 

In order to demonstrate the involved physics, we shall first 
derive the dispersion relation for surface plasmon-polaritons 
ifTOll in a simplified structure. Let us consider an electro- 
magneic wave in the vicinity of a plane interface between 
two semi-infinite media, a metal and a QD/dielectric compos- 
ite, and assume that the wave is p— polarized. We choose x 
axis along the direction of propagation of the electromagnetic 
wave within the interface plane and z axis perpendicular to the 
interface. Assuming the electromagnetic field's time- and x- 
coordinate dependence to be ~ exp {ikx — iujt) (where k and 
w are the wavenumber and frequency), we can write down the 
wave equation in the form: 
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It possesses solution corresponding to an evanescent wave, 

ii;«(z)=£;W(0)exp(-piz), (9) 
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where p,„ ~ \/k^ — uj'^em/c^', m = 1,2 and Rc(p,„) > 
as it is characteristic of surface waves. Let to = 1 stand 
for the metal and ei = eni'^) is the Drude-type dielectric 
function, and to = 2 corresponds to the composite medium, 
£2 = ^*{^)- By applying boundary conditions at the interface 
we obtain the SPP dispersion relation in the form: 

P2 Pi 

It can be solved for the wavenumber k, yielding 
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FIG. 2: (color on line) Panels (a, c, e): Real (solid lines) and imag- 
inary (dashed lines) parts of the SPP wavevector versus frequency, 
calculated for the interface between semiinfinite CdSe QD - PMMA 
composite and silver; Panels (b, d, f): reflectivity R versus angle 
of incidence 9, and frequency cj for the ATR structure with a glass 
prism (Sg — 2.9584), a silver film of thickness d = 53.3 nm, and 
a semiinfinite composite medium. The composite dielectric function 
was calculated using MMGA with a = 3 nm and Aa = (panels a, 
b), Aa = 0.005a (panels c, d), and Aa = 0.1a (panels e, f). 



The frequency dependence of the real and imaginary parts of 
k for an Au/QD-PMMA interface is shown in Figs. |2|a), 
(c), and (e). Im(fc) increases drastically in the vicinity of 
resonance frequencies that are determined by the relation 
Re {e*{uj) + eM(w)) — and are quite close to the QD ex- 
citon transition frequencies. It reflects the resonant coupling 
between SPPs and QD excitons. Again, the resonance peak 
broadens with the increase of the QD size dispersion because 
of the weaker coupling involving fewer QDs for each given u. 

As known, SPPs can be probed in the ATR geometry 
schematically shown in in Fig. |3] where the frquency and the 
x-projection of the wavevector of an external electromagnetic 
wave can match uj and Re (A:) of the surface polaritons 1.10,1 . 

sin 9 , 

where eg is the dielectric constant of the prism and 6 is tha 
angle of incidence. The matching is achieved by adjusting 
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FIG. 3: ATR structure composed of a glass prism, a silver plate, 
and two composite layers containing semiconductor NCs of different 
average size 



u! and/or 9 and is detected by measuring the reflectance, R, 
which shows characteristic dips corresponding to the reso- 
nance transfer of the electromgnetic energy into surface po- 
laritons. We have performed calculations of R for the ATR 
structure consisting of a prism, an Ag film |13|, and one or 
two layers of CdSe QD-PMMA composites (Fig. |3]l, which 
are similar to those outlined above for the simplified struc- 
ture containing just two semi-infinite media. The minima in 
the ATR spectra can be observed in Figs. |2|b),(d), and (f). 
They resemble the structure of the corresponding Re(fc) vs u 
dependencies [Figs. [2]; a), (c), and (e)]. 



RESULTS AND DISCUSSION 

Let us consider the structure of the SPP-exciton resonance 
in more detail. In the case of semiinfinite composite with 
CdSe NCs [Figs. [4]^a)-(c)] the resonant interaction of modes 
occurs at O w 68° for the chosen set of parameters. In the 
vicinity of this angle the dependence exhibits two min- 
ima with approximately equal depth, while for other angles 
of incidence one of the minima is significantly deeper than 
the other [see Fig. |4|b)]. This is characteristic of mode anti- 
crossing and corresponds to the experimental observation of 
Ref 4 The positions of these minima (i.e., the angles of in- 
cidence, Omin, at which the reflectance reaches its minimum, 
Rmin, for a given frequency) are depicted in Fig. |2|c), where 
one can clearly see the resonant SPP-exciton interaction, or 
mode anticrossing. At the same time the depth of the reso- 
nant minima of Rmin{^) at the corresponding angles dmin 
demonstrates that there are at least three points {oj « 2.36 eV, 
w « 2.6 eV, w « 2.72 eV) where the reflectivity of the ATR 
structure reaches zero, Rmin = 0. 

How does the CdSe QD-PMMA layer thickness influence 
the ATR spectrum? First, the ATR resonance minima are 
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FIG. 4: (color on line) Reflectivity versus angle of incidence 6 and 
frequency uj [panels (a, d, g)]; reflectivity versus frequency uj for a 
fixed 9 [panels (b, e, h)]; minimal reflectivity Rmin (dashed lines) 
and corresponding angle of incidence 9min (solid lines) for a given uj 
[panels (c, f, i)] for the ATR structure with glass prism Eg = 2.9584, 
silver film of thickness d = 53.3 nm, and CdSe QD-PMMA com- 
posite layer of different thickness: semiinfinite (a-c), di — 53.3 nm 
(d-f) and di = 19.78 nm (g-i). In the case of finite di, the medium 
at z < —di is vacuum. The composite dielectric function was calcu- 
lated using MMGA with / = 0.1, a = 3 nm and = 0.1a. 

shifted to lower angles of incidence [compare Figs. Qa) and 
(d)]. Secondly, the anti crossing of modes becomes weaker as 
manifested by the smaller separation of the reflectance min- 
ima in Fig. |4|e) for 9 w 50°, as well as in Fig. |4|f)]. 
Thirdly, the values of Rmin far from the exciton resonance 
frequency are higher than in the case of semi-infinite compos- 
ite medium. Further decrease of the CdSe QD-PMMA layer 
thickness [Figs. |4|g)-(i)] leads to the damping of the SPP- 




FIG. 5: (color on line) Reflectivity (a, d), and relative electric field 
square amplitude \Ex{z)/ Ei\'^ at z = — d (b, e) and z — —{d + di) 
(c, f) versus angle of incidence and frequency for the ATR struc- 
ture consisting of glass prism Eg = 2.9584, Ag film of thickness 
d = 53.3 nm and two CdSe QD-PMMA composite layers, one of 
thickness di = 19.78 nm (a-c) or di — 53.3 nm (d-e) and average 
QD radius ai = 3.7 nm, and the other semi-infinite, with average 
QD radius a2 — 3 nm. The dielectric functions for both composite 
layers were calculated using MMGA with / = 0.1 and QD radius 
dispersion Aa = 0.1a in all cases. 

exciton interaction. For example, in the case of Fig. |4|g) the 
resonance is hardly distinguishable (it is confirmed also by 
Fig. |4ji) where the positions of resonant minima are shown 
explicitly). Also from Fig. |4|h) one can see that for the angles 
6 w 41° and 9 w 41.5° two reflectivity minima practically 
merge and cannot be resolved. 

Now let us turn to the structures with the same silver film 
and two composite layers, one of thickness c?i containing QDs 
of average radius ai and the other (semi-infinite) with QDs of 
radius 02 (Fig. [Slwith d2 00). In the ATR spectrum of this 
structure [Figs. |5ja), (d)], one can observe resonances corre- 
sponding to both ai [hu! w 2.2 eV and hu> 2.7 eW in Figs. 
0:a), (d)] and 03 QDs [huj « 2.5 eV in Figs. [Sfa), (d)]. In- 
creasing di results in a stronger anticrossing between the SPP 
and the ai QD exciton modes, and in a weakeing of the 02 
resonance [compare Figs. |5ja) and (d)]. Nevertheless, choos- 
ing an appropriate composite layer thickness, one can achieve 
SPP coupling to QD excitons localized in the composite layer 
which is not adjacent to the metallic film. 

ATR minima correspond to some particular values of u and 
9 for which the energy of the incident electromagnetic wave 
is most efficiently transferred to SPP and, consequently, to the 
QD excitons. This is shown in Figs. |5]^b),(c),(e) and (f), where 



5 




■ CO, eV 



(s, eV' 



FIG. 6: (color on line) Reflectivity versus angle of incidence and fre- 
quency for the ATR structure consisting of glass prism Eg — 2.9584, 
Au film of thickness d = 53.3 nm and two CdSe QD-PMMA 
composite layers, one of thickness di and average QD radius ai 
barai = 4.2876 nm (c, d) and the other semiinfinite with average 
QD radius d2 = 3 nm. The following parameters have been used: 
di = 53.3 nm, oi — 3nm 02 — 3.2 nm (a); di = 79.12 nm, 
ai — 3nm 02 = 3.2 nm (b); di — 19.78 nm, ai — 4.2876 nm 
a2 = 3nm (c); di = 53.3 nm, ai = 4.2876 nm a2 = 3nm (d). 
Insets in panels (a, b) present the dependence of the reflectivity on 
frequency for a fixed angle of incidence. 



the squared average amplitude of the electric field is plotted 
for each of the CdSe QD-PMMA layers. It represents the 
intensity of excitation of QD luminescence, whose character- 
istic colour is determined by the QD size. Therefore one can 
adjust the relative PL intensity of different colours emitting by 
two or more composite layers containing QDs of different av- 
erage radii. It opens the possibility of building a lightning de- 
vice with colour characteristics which can be controlled either 
by the excitation frequency or by the ATR incidence angle. 

The average sizes of the NCs embedded in different lay- 
ers can be chosen in such a way that the frequencies of the 
excitonic resonances are close to each other. The ATR spec- 
tra calculated for this situation are depicted in Figs. |6]^a),(b). 
As it is can be seen from the insets, in this case (that can be 
called double resonance) it is possible to achieve the situation 
where three reflectivity minima coiTespond to one value of 6. 
We also checked another special situation where the higher 
energy exciton transition in ai QDs matches the ground state 
transition in the 02 dots [see Figs. |5|c), (d)] but apparently no 
new effects arise in this case and there are only two resonance 
bands (at huj w 2.1 eV and huj ^ 2.5 eV) seen. 

Finally, we considered still another potentially interesting 
sandwich-type structure containing two [Fig. [Tj^a)] or three 
[Fig. |7];b)] interfaces between silver and CdSe QD-PMMA 
composite. For such structures, the calculated reflectivity 
spectra contain two or three minima (for a gived frequency), 
respectively. One can observe a rather sophisticated "bunch- 
ing" of modes in the vicinity of the excitonic resonance. 

In conclusion, we have shown that the resonance cou- 
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FIG. 7: (color on line) (Panel a) Reflectivity versus angle of inci- 
dence and frequency for the ATR structure consisting of the glass 
prism, one silver film of thickness d — 53.3 nm, one CdSe QD- 
PMMA composite layer with average QD radius ai — 3nm and 
thickness di = 300 nm, and a semiinfinite silver layer (a) or a sec- 
ond silver film of thickness d\^ — 69.23 nm and a semi-infinite 
CdSe QD-PMMA composite with average QD radius ai (b). 



pling between the plasmon-polaritons propagating along the 
metal/NC-layer interface and excitons confined in chemically 
synthesized semiconductor NCs, experimentally observed in 
Ref. 0] can be described theoretically using the appropriate 
effective dielectric function for the NC composite layer and 
standard multilayer optics. The SPP-exciton interaction pro- 
duces a considerable effect on the optical properties of the 
structure if the dispersion of the NC size in the composite 
layer is not too large. We have shown that combining sev- 
eral composite layers with appropriately sized quantum dots 
and/or more than one metallic films can result in interesting 
interactions between the various SPP and exciton modes. Ow- 
ing to these interactions, the energy of an incident electromag- 
netic wave can be distributed, by means of surface plasmons, 
between the different QD species, as it has been suggested for 
molecules adsorbed on a metallic surface |[T2l . In particular, 
it opens the possibility to control the relative intensity of light 
of different colors, emitted by the QDs of different sizes. The 
work was supported by the FCT (Portugal) under the grant 
PTDC/FIS/113199/2009. 
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